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Abstract 

Modelling  and  hardware-in-the-loop  (HIL)  testing  of  fuel  cell  components  and  entire  systems  open  new  ways  for  the  design  and  advance 
development  of  FCs.  In  this  work  proton  exchange  membrane  fuel  cells  (PEMFC)  are  dynamically  modelled  within  MATLAB-Simulink  at  various 
operation  conditions  in  order  to  establish  a  comprehensive  description  of  their  dynamic  behaviour  as  well  as  to  explore  the  modelling  facility  as  a 
diagnostic  tool.  Set-up  of  a  hardware-in-the-loop  (HIL)  system  enables  real  time  interaction  between  the  selected  hardware  and  the  model. 

The  transport  of  hydrogen,  nitrogen,  oxygen,  water  vapour  and  liquid  water  in  the  gas  diffusion  and  catalyst  layers  of  the  stack  are  incorporated 
into  the  model  according  to  their  physical  and  electrochemical  characteristics.  Other  processes  investigated  include,  e.g.,  the  membrane  resistance 
as  a  function  of  the  water  content  during  fast  load  changes.  Cells  are  modelled  three-dimensionally  and  dynamically.  In  case  of  system  simulations 
a  one-dimensional  model  is  preferred  to  reduce  computation  time.  The  model  has  been  verified  by  experiments  with  a  water-cooled  stack. 

©  2006  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

By  means  of  mathematical  modelling  the  development  and 
design  of  fuel  cell  systems  can  be  significantly  accelerated. 
To  optimise  and  develop  the  fuel  cell  hardware  first  a  detailed 
description  of  a  PEMFC  three-dimensional  model  is  given.  In 
case  of  system  simulations  and  HIL  applications  with  regard  to 
computational  time  a  one-dimensional  PEMFC  model  is  used. 
Both  the  three-dimensional  as  well  as  one-dimensional  model 
are  implemented  in  MATLAB/Simulink. 

The  main  model  outputs  are  cell  voltage,  temperature  as 
well  as  electrical  and  thermal  power,  which  are  calculated  with 
respect  to  different  operating  conditions,  such  as  load  cur¬ 
rent,  pressure,  flow,  gas  composition,  moisture,  etc.  The  model 
uses  basic  electrochemical  and  physical  principles  with  special 
emphasis  on  the  description  of  dynamic  behaviour. 

The  vehicle  model  consists  of  the  commercially  available 
simulation  tool  CarSim  where  the  internal  combustion  engine  is 
replaced  by  a  model  of  an  asynchronous  motor  and  an  inverter. 
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After  the  model-based  analysis  the  hardware-in-the-loop 
(HIL)  vehicle  system  is  introduced,  which  allows  to  investi¬ 
gate  the  behaviour  of  the  real  components  in  a  virtual  envi¬ 
ronment.  The  main  components  of  the  fuel  cell  HIL  system, 
their  properties  will  be  described  in  detail  below,  are  the  com¬ 
mand  station,  target  node  PC  with  real  time  operating  system, 
signal-conditioning  unit,  device  under  test  and,  for  presentation 
purposes,  a  separate  human  machine  interface  consisting  of  a 
steering-wheel  and  a  break  and  throttle  pedal. 

2.  Three-dimensional  fuel  cell  model 

The  dynamic  PEFC  three-dimensional  model  (Fig.  1) 
includes  the  flowfield,  gas  diffusion  layer  (GDL),  catalyst  layer 
of  the  anode  and  cathode  as  well  as  the  membrane.  The  model 
has  been  implemented  in  MATLAB/Simulink  software  and  the 
simulation  can  run  with  discrete  or  continuous  time  steps. 

The  number  of  elements  in  each  direction  can  be  chosen 
freely  with  exception  of  the  catalyst,  which  is  modelled  as 
a  border  layer  and  the  flowfield.  The  flowfield  is  modelled 
only  in  x-y  direction,  because  the  gases  have  to  follow  the  gas 
channels  and  so  the  gas  composition  is  changing  mainly  along 
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Nomenclature 

^An 

Anode  charge  transfer  coefficient  (-) 

Q'Ca 

Cathode  charge  transfer  coefficient  (-) 

£gdl 

Porosity  GDL  (-) 

M 

Potential  (V) 

X 

Membrane  water  content  (-) 

A 

Cell  area  (m2) 

G 

Molar  concentration  of  component  i  (molm-3) 

d 

Width  of  the  membrane  (m) 

Di 

Diffusive  coefficient  of  component  i  (-) 

E 

Activity  coefficient  (-) 

F 

Faraday  constant  (C  mol-1) 

A  G° 

Gibbs  energy  (J  mol-1) 

i 

Current  density  (A  cm-2) 

70,  An 

Exchange  current  density  anode  (A  cm-2) 

jo,  Ca 

Exchange  current  density  cathode  (A  cm-2) 

Ji 

Diffusive  flux  of  component  i  (molm-2  s-1) 

k 

Proton  conductivity  (S  m-1) 

Mi 

Molar  mass  (kg) 

ni 

Mass  (mol) 

m 

Mass  flow  rate  (mol  s-1) 

Ni 

Molar  flux  of  component  i  (mol  m-2  s-1) 

Pi 

Partial  pressure  of  component  i  (pa) 

R 

Gas  constant  (j  mol-1  K-1) 

7?ei 

Electric  resistance  (£2  cm2) 

^mem 

Ionic  resistance  (£2  cm2) 

Swelling  factor  (-) 

t 

Time  (s) 

t\v 

Transport  number  (-) 

T 

Cell  temperature  (K) 

T° 

Reference  temperature  (K) 

Uo 

Open  cell  voltage  (V) 

uAn 

Anode  voltage  losses  (V) 

UC a 

Cathode  voltage  losses  (V) 

f^cell 

Cell  voltage  (V) 

Uq 

Ohmic  voltage  loss  (V) 

um 

Mixture  velocity  (ms-1) 

V 

Volume  (m3) 

X 

Mol  fraction  (-) 

z 

Length  (m) 

this  direction.  The  flow-field  could  be  a  linear  one  or  a  meander 
type  and  it  is  possible  to  chose  the  direction  of  the  anode  and 
cathode  flow  in  linear,  counter  or  cross.  The  pressure  and  the 
gas  composition  is  changing  along  the  channel  depending  on 
the  current  density  in  each  simulated  element. 

The  gas  diffusion  layer  incorporates  the  gas  transport  from 
the  flowfield  to  the  catalyst  layer  and  the  removal  of  the  reactants. 
Here  the  transport  of  liquid  water  is  also  included.  For  the  cat¬ 
alyst  layer  the  main  electrochemical  processes  are  considered. 
The  mentioned  fuel  cell  components  and  especially  the  mem¬ 
brane  and  their  influence  on  the  performance  could  be  analysed 
in  detail. 


Fig.  1.  Fuel  cell  model. 


2.7.  Flowfield 

After  inlet  the  gases  have  to  follow  the  channels  through  the 
flowfield  plate.  For  computational  convenience  the  gases  are 
supposed  to  obey  the  ideal  gas  law.  At  the  anode  H2,  H2O  and  at 
the  cathode  O2,  N2,  H2O  are  considered.  The  pressure  distribu¬ 
tion  of  all  gas  components  and  the  whole  pressure  is  calculated 
along  the  gas  channels.  This  is  influenced  by  the  transported 
species  via  concentration  and  pressure  gradients  of  the  processes 
in  the  GDL,  which  are  explained  in  detail  in  the  next  paragraph. 

2.2.  Gas  diffusion  layer 

Gas  diffusion  layers  serve  as  current  collectors  allowing  ready 
access  of  the  fuel  and  oxidant  to  anode  and  cathode  catalyst 
surfaces,  respectively.  A  gas  diffusion  layer  usually  contains  a 
mixture  of  carbon  black  powder,  hydrophobic  dispersion  agent, 
and  solvents,  applied  onto  carbon  paper  or  woven  carbon  cloth 
substrates  (both  have  macropores),  to  form  a  uniform  thin  micro- 
porous  layer  on  top  of  the  macroporous  layer  of  the  substrates. 
Here,  gas  diffusion  layers  are  manufactured  using  Toray  papers 
as  the  macroporous  substrate. 

In  the  model,  the  main  input  parameters  of  the  GDL  anode 
and  the  GDL  cathode  are: 

•  i,  current  density, 

•  T,  temperature, 

•  Pft2,  Po2,  Pn2>  Ph2o,  partial  pressures,  and 

•  geometric  and  material  characteristics. 

In  the  following,  the  calculations  of  the  GDL  on  the  basis  of 
reference  [1]  will  be  explained.  The  diffusive  transport  is  mod¬ 
elled  using  Fick’s  first  law.  However,  it  is  strictly  valid  only  for 
binary  mixtures  and  for  some  mixtures  of  three  or  more  com¬ 
ponents,  the  Maxwell-Stefan  equations  are  required  to  obtain 
reliable  results. 

The  governing  equations  for  the  gas  mixture  will  be  described 
by  the  concentration  q  in  the  subsequent  equations.  The  next 
table  presents  the  values  that  q  takes  in  the  anode  and  in  the 
cathode.  The  conservation  of  mass  for  the  gas  mixture  takes 
the  form  of  the  Eq.  (1).  The  evaluation  is  made  using  the  ODE 
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Matlab  functions. 


dCi 

~dt 


=  -VNi  =  -V(Hc/) 


(1) 


where  um  is  the  molar- averaged  mixture  velocity  and  the  diffu¬ 
sive  flux  (measured  relative  to  the  molar- averaged  velocity).  Eq. 
(1)  describes  the  total  convective  plus  diffusive  flux  of  each  gas. 
The  diffusive  flux  is  given  by  Fick’s  law,  which  states  that  the 
flux  of  one  component  relative  to  the  molar-averaged  velocity  is 
proportional  to  the  gradient  in  mole  fractions  via, 


Ji  =  A  CiDiA 


£gdl 

Az 


(2) 


where  the  diffusion  coefficient  D[  is  given  by, 


Di  = 


1  - 


1 


PMIX  )  y^pk/ ffMIx) 


(3) 


k 


The  second  term  of  Eq.  (1)  is  the  convective  term: 

£gdl 


CiUm  =  AcjApMixA 


A  z 


(4) 


In  both  Eqs.  (2)  and  (4)  the  porosity  £gdl  is  represented  by: 
3TH20^H20lici 


£GDL  = 


(5) 


Mh2o  and  p±i2Q\iq  are  the  molar  mass  and  the  density  of  the 
water  respectively,  and  finally  is  the  A  V  volume  of  the  GDL.  The 
pressure  of  the  mixture  pmix  is  evaluated  assuming  the  gas  to 
be  ideal  and  adiabatic.  It  depends  linearly  on  the  concentration. 


Pmix  =  Y ^Pi 


(6) 


where  pi  is  the  partial  pressure  of  each  gas, 
Pi  =  ctRT 


(7) 


R  is  the  universal  gas  constant  and  T  the  temperature.  Although 
there  are  significant  temperature  variations  within  a  fuel  cell,  the 
model  considers  the  system  to  be  isothermal. 


^h2o  =  —  A  (10) 

zr 

where  i  is  the  current  density  and  A  the  area  of  the  membrane. 
2.4.  Membrane 

The  amount  of  water  fed  into  the  cell  by  humidified  reactant 
gases  depends  on  the  gas  temperature,  pressure  and  flow  rate, 
whereas  water  production  on  the  cathode  depends  mainly  on  the 
current  density,  to  which  it  is  directly  proportional.  In  a  PEM  fuel 
cell  water  is  transported  along  the  gas  channels,  but  also  through 
the  membrane  and  the  electrodes.  Transport  processes,  which 
are  schematically  represented  in  Fig.  2,  occur  in  an  ionomere 
membrane. 

For  the  water  transport,  the  principles  driving  modelled  forces 
are: 

•  Convective  force 

•  Osmotic  force  or  diffusion 

•  Electric  force 


The  first  of  these  forces  arises  from  a  pressure  gradient,  the 
second  from  a  concentration  gradient,  and  the  third  from  the 
migration  of  protons  from  anode  to  cathode  and  their  effect 
(drag)  on  the  dipole  water  molecules.  For  water  and  protons 
the  mass  conservation  equation  can  be  represented  as: 

-f  =  -VNi,  i  =  H20,  H+  (11) 

at 

where  c/  is  the  molar  concentration  and  A/  the  molar  flux  due  to 
electro- osmotic  driving  forces  and  convection.  The  molar  flux 
can  be  represented  as: 

Ni  =  Ji  +  Cium  (12) 


where  um  is  the  mixture  velocity  and  Ji  the  diffusive  flux.  In  the 
model,  the  diffusive  flux  for  the  water  and  the  proton  current  is 
calculated  on  the  bases  of  [2]  as: 


twk 

^2  dz 


DujCH20  9mh2o\  1 
RT  dz~ )  1+sX 


(13) 


2.3.  Catalyst  layer 


For  the  catalyst  layer  the  main  gas  transport  equations  as 
in  the  GDF  are  included.  A  sink  term  is  used  to  account  for 
depletion  of  oxygen  at  the  cathode  and  depletion  of  hydrogen 
at  the  anode.  The  electrochemical  and  physical  processes  are 
complex  at  the  catalyst  layer.  First-order  kinetics  is  assumed 
using  the  Bulter-Volmer  equation.  For  the  liquid  water  mass 
conservation  equation,  the  source  terms  at  the  cathode  and  anode 
are  respectively: 


hiin  =  — A 
2  2  F 


i 

4 F 


A 


Anode  GDL  Membrane  Cathode  GDL 


V 

IF 

1  : 

A1  \ 

/  Convective  force,  HzO 
Osmotic  force,  H20 

\J  1/ 

Electric  force,  H20,  H  ' 

T  T  i 

Anode  catalyst  Cathode  catalyst 

Fig.  2.  Direction  of  the  water  transport  forces  inside  the  membrane. 
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■Ar  — 


k  9/Xj^ 


+ 


twk  dji U2o\  1 

)  l +sX 


F 2  dz  F2  3  z 

The  potential  gradient  from  water  is  calculated  through: 

3mh2o  _  RT  3xh2o 
XR20  dz 

Inserting  the  Eqs.  (14)  and  (15)  in  (13) 

(  t2k  RT  3xh9o 

7h2o  =  (  —  +  sF)  +  ~Fi - A — 

V  F  *h2o  dz 

_  DwCh2o  dxu2o\  1 

*H20  dz  J  1  +  sk 

where  the  transport  number  for  Nation®  is  evaluated  as: 
A 


(14) 


(15) 


(16) 


tw  =  2.5  — 

22 

(17) 

where  in 

,  _  «H2o 

(18) 

^so3h 

withnso3H  =  1.8634  x  103  mol/m2 

The  swelling  factor  takes  the  value  s  =  0.0123. 

The  temperature  dependence  of  k  is  calculated  with: 

k  =  j^e-EAMim-a/T^) 

where  Tre f  =  353.15  K  and  Kre f  is  given  by: 

£ref  =  0.0013A3  +  0.0298A2  +  0.2568A 


(19) 


(20) 


In  (19)  the  coefficient  of  activity  appears,  which  is  calculated 
as: 

£a,k  =  2640e“oa  +  1183  (21) 

The  diffusion  coefficient  is  calculated  by  the  expression: 

Dw  =  2Ue(_2436/r)  (22) 

The  molar  fractions  are: 

A 


xu2o  = 


xn+  — 


A  +  2 

1 


(23) 


A  +  2 

The  concentration  of  the  water  is: 
^h2o 


ch2o  = 


V 


(24) 


where  V  is  the  volume  of  the  membrane.  The  current  is  given  by 
the  flow  of  protons: 


2.5.  Cell  voltage 


The  cell  voltage  is  described  by  the  following  equation: 


UcqM  =  Uo  —  U An  —  Uq  —  Uqsl 


(27) 


where  Uo  is  the  open  cell  voltage,  U\n  and  Uca  are  the  voltage 
losses  on  anode  and  cathode  side,  respectively,  and  Uq  repre¬ 
sents  the  ohmic  losses  due  to  ionic  conductivity  of  the  membrane 
and  due  to  electronic  conductivity  of  bipolar  plates,  diffusion  and 
catalyst  layers. 

The  open  cell  voltage  is  a  function  of  temperature  and  partial 
pressures  of  reactants  and  is  calculated  as  follows: 


Uo  =  ~ 


A  G°  A  S 
+  —  (T 


-  T°)  + 


2  F 


IF 


RT 

~2F 


1  ,  PO2  .  -  Pn2 




F 


F 


(28) 


Based  on  the  Bulter-Volmer  equation  [3],  which  describes  the 
reduction  and  oxidation  processes  on  the  electrode,  and  neglect¬ 
ing  the  opposite  reaction  on  each  electrode,  the  losses  through 
the  charge  transfer  can  be  described  as  follows: 


RT  i 

(O' An  +  OiQ a)F  z'o,An 


(29) 


RT 

OiCaF 


In 


i 

*0,Ca 


(30) 


The  exchange  current  density  io  and  transfer  coefficient  a  are 
measures  for  the  reaction  kinetics.  The  exchange  current  density 
is  concentration  and  temperature  dependent.  Since  the  exchange 
current  density  for  the  hydrogen  electrode  is  considerably  higher 
than  oxygen,  charge  transfer  losses  at  the  anode  are  small  and 
may  be  neglected. 

Ohmic  losses,  mainly  associated  with  the  membrane  part,  are 
given  by  the  following  equation: 


Uq  =  Re\1  +  RMemi 


(31) 


where  Rei  represents  the  losses  through  the  electronic  conduc¬ 
tivity  and  Ryiem  the  voltage  loss  due  to  ionic  conductivity  of 
the  membrane.  The  conductivity  of  the  membrane  depends  on 
temperature  and  water  content  [4] . 

Because  of  the  low  charge  transfer  resistance  the  resulting 
time  constant  is  also  small.  Therefore,  the  double  layer  capaci¬ 
tance  on  the  anode  side  can  be  neglected. 


2.6.  Simulation  results 


■V  =  f  (25) 

F 

The  resistance  of  the  membrane  corresponds  to: 
d 

Ryiem  —  T  (26) 

k 

where  d  is  the  width  of  the  membrane  and  k  is  calculated  by  Eq. 
(19). 


The  dynamic  results  of  the  3D  fuel  cell  model  are  shown  in 
this  chapter  and  compared  with  the  experimental  results.  The 
measurement  was  done  with  a  water-cooled  stack  within  five 
cells.  The  active  area  was  100  cm2  and  the  membrane  type  used 
was  GORE  PRIMEA  MESGA  with  a  thickness  of  35  p,m.  The 
GDL  employed  was  SIGRACET  GDL,  GDL-10  BB,  420  pan. 
The  meander  type  flowfield  with  a  depth  of  1 .3  mm  was  used  in 
parallel  counter  flow. 
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Fig.  3.  Partial  pressure  of  oxygen  and  water  in  the  flowfield. 


Fig.  5.  Membrane  water  content. 


2.6.1.  Influence  on  partial  pressure 

The  partial  pressure  of  oxygen  and  the  water  vapour  pressure 
in  the  flowfield  are  shown  in  Fig.  3  for  a  current  density  of 
0.8  A  cm-2.  The  calculation  was  done  with  10  elements  in  v  and 
y  direction.  In  this  case  the  oxygen  pressure  is  decreasing  along 
the  gas  channel  of  the  meander  from  the  inlet  to  the  outlet.  The 
water  vapour  pressure  is  increasing  from  the  inlet  to  the  outlet. 

The  GDL  is  simulated  within  several  layers.  As  shown  in 
Fig.  4  the  oxygen  partial  pressure  in  the  layer  next  to  the  cata¬ 
lyst  layer  is  always  higher  than  0.1  bar  at  a  current  density  of 
0.8  A  cm-2.  The  oxygen  pressure  is  increasing  according  to  the 
flowfield.  During  these  conditions  there  was  also  liquid  water  in 
the  GDL.  However,  the  amount  was  low  and  did  not  lead  to  a 
significant  reduction  of  reactive  gases. 

2.6.2.  Membrane  performance 

The  performance  of  the  cell  is  strongly  influenced  by  the 
resistance  and  the  water  content  of  the  membrane.  The  equa¬ 
tions  of  the  membrane  model  are  based  on  Nation®  as  well 
known  material.  The  parameters  of  the  used  Gore  membrane 


are  different  to  Nation®  and  so  during  the  calculations  a  param¬ 
eter  fit  was  necessary. 

In  Fig.  5  the  amount  of  water  is  presented  for  three  different 
times  during  the  simulation  in  the  x-z  layer.  The  amount  of  water 
before  the  current  step  at  0.05  A  cm-2  is  almost  equal  in  the 
membrane  and  low,  because  of  minor  water  production.  Twenty 
seconds  after  a  current  step  the  amount  of  water  is  increasing 
because  of  the  high  current  density  (0.80  A  cm-2).  After  50  s, 
the  transients  are  almost  decayed  and  the  water  content  of  the 
membrane  is  higher  near  the  cathode  than  at  the  anode. 

The  change  of  the  membrane  resistance  in  x—y  direction  for 
the  same  current  step  is  shown  in  Fig.  6.  The  influence  of  the 
meander  flow-field  on  the  water  amount  leads  to  high  differences 
in  the  membrane  resistance.  The  average  value  of  the  resistance 
before  the  current  step  is  about  five-times  higher  than  50  s  after 
the  current  step. 

The  current  density  is  not  constant  because  of  an  inhomoge¬ 
neous  resistance  in  x-y  direction  (Fig.  6).  At  an  average  density 
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Fig.  4.  Partial  pressure  of  oxygen  and  water  in  GDL. 


Fig.  6.  Resistance  of  the  membrane. 
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of  0.05  A  cm-2  the  current  density  (Fig.  7)  is  not  varying  a  lot. 
At  higher  currents  the  membrane  acquires  a  high  local  current 
density.  At  0.80  A  cm-2  the  highest  local  current  density  is  up 
to  1.00  A  cm-2  and  the  lowest  is  only  about  0.50  A  cm-2.  This 
leads  to  higher  voltage  losses  at  the  membrane  and  a  lower  cell 
voltage.  This  influence  is  dynamically  described  in  the  next  para¬ 
graph. 

2.6.3.  Dynamic  behaviour 

The  comparison  of  the  dynamic  behaviour  of  the  fuel  cell 
model  and  the  measurement  are  presented  for  a  relative  humidity 
at  the  cathode  of  80%  and  at  the  anode  of  20%.  The  stoichiom¬ 
etry  was  4  at  the  cathode  and  2  at  the  anode.  The  cell  voltage  of 
the  simulation  is  compared  in  Fig.  8  with  the  measurement  for 
dynamic  loads.  The  presented  voltages  are  standardised.  The  fig¬ 
ure  shows  always  an  excellent  agreement  between  measurement 
and  simulation. 

After  a  current  step  to  higher  values  the  voltage  is  at  first 
decreasing  sharply,  then  increasing  in  the  next  60  s  to  its  stable 
value.  This  behaviour  can  be  explained  by  the  humidity  change 
of  the  membrane,  which  is  increasing,  because  of  the  water  pro¬ 
duction  during  high  current.  Switching  back  to  lower  current 
values  the  membrane  is  drying  out  and  the  voltage  is  decreasing 
after  the  first  few  seconds. 


Time  in  Min 

Fig.  8.  Comparison  measurement  and  simulation. 

As  shown  in  Fig.  9  the  whole  fuel  cell  system  can  only  be 
tested  if  all  sub-systems  and  units  are  available.  In  contrast, 
Fig.  10  shows  the  configuration  if  there  is  a  possibility  to  inte¬ 
grate  HIL  simulation. 

As  demonstrated  in  Fig.  10  it  is  now  possible  to  replace  defec¬ 
tive  or  unavailable  components  with  their  mathematical  models 
and  perform  analysis  of  operational  behaviour  of  the  fuel  cell 
system. 

Because  of  the  possibility  to  operate  the  hardware  together 
with  simulation  models  one  of  the  main  aspects  for  the  reali¬ 
sation  of  HIL  test  stands  are  the  safety  requirements.  A  high 
safety  standard  is  reached  through  the  implementation  of  mul¬ 
tistage  switch  off  procedure  on  the  PLC  (Fig.  10).  Furthermore, 
a  plausibility  check  is  done  by  the  PLC  for  the  set  values  of  the 
controllers  calculated  in  the  model,  before  being  passed  on  to 
the  hardware.  The  PLC  will  restrict  the  set  values  to  reasonable 
limits  and  prevent  dangerous  operation  conditions.  This  is  espe¬ 
cially  important  in  case  of,  e.g.,  mixing  hydrogen  and  oxygen 
on  the  anode  side,  or  when  calculating  the  pressure  set  values 
by  the  model. 


3.  HIL  test  benches 

For  the  development  of  fuel  cell  stacks  or  fuel  cell  systems, 
there  is  an  increasing  need  in  suitable  test  benches.  The  MAG¬ 
NUM  company  has  been  one  of  the  first  suppliers  of  fuel  cell 
test  benches  and  in  particular  of  HIL  capable  test  benches.  The 
advantage  of  using  HIL  test  stands  is  that  already  in  an  early  stage 
of  development  testing  of  whole  system  can  be  performed,  even 
though  not  all  the  hardware  is  available.  The  following  Fig.  9 
shows  the  configuration  of  a  conventional  test  stand. 


Fig.  9.  Conventional  test  bench. 
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Fig.  10.  HIL  capable  test  bench. 


Another  important  aspect  regarding  automotive  applications 
is  the  dynamic  operation  of  the  test  bench.  Therefore,  fast  mass 
flow  and  pressure  controllers  are  integrated  and  also  fast  data 
acquisition  and  data  exchange  between  the  PLC  and  control  PC, 
respectively,  simulation  PC  are  realised. 

In  order  to  afford  dynamic  changes  of  the  gas  humidification 
in  the  test  stand,  model  based  control  algorithms  are  imple¬ 
mented.  Thermodynamic  simulations  have  been  integrated  to 
develop  new  control  algorithms  for  temperature  and  moisture. 
These  methods  allow  to  determine  the  actual  moisture  with  high 
accuracy,  by  measuring  solely  temperature  and  considering  heat 
capacities  of  gases  and  equipment. 

Fig.  1 1  gives  an  overview  of  the  general  configuration  of  the 
above  described  test  bench  and  the  interaction  of  the  test  stand 
with  the  control  and  simulation  PC. 

As  shown  in  Fig.  1 0  the  fast  data  transfer  between  the  automa¬ 
tion  PC,  HIL  computer  and  the  PLC  is  realised  via  ethernet 
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Fig.  11.  Automation  and  data  transfer  within  HIL  test  bench. 


connection.  The  HIL  simulation  can  be  started  by  the  operator 
choosing  an  adequate  operation  mode  on  the  automation  PC.  The 
simulation  model  is  running  on  the  HIL  PC  in  real  time  writing 
selected  set  values  to  the  PLC  and  simultaneously  reading  actual 
values  from  PLC. 

3.1.  HIL  fuel  cell  vehicle  simulation  system 

In  case  of  system  simulations  and  HIL  applications  with 
regard  to  computational  time  a  one-dimensional  PEMFC  model 
is  used,  which  is  implemented  in  MATLAB/Simulink.  The  main 
model  outputs  are  cell  voltage,  temperature  as  well  as  electrical 
and  thermal  power,  which  are  calculated  with  respect  to  differ¬ 
ent  operating  conditions,  such  as  load  current,  pressure,  flow, 
gas  composition,  moisture,  etc.  The  vehicle  model  consists  of 
the  commercially  available  simulation  tool  CarSim  where  the 
internal  combustion  engine  is  replaced  by  a  model  of  an  asyn¬ 
chronous  motor  and  an  inverter. 

After  the  model  based  analysis  the  HIL  vehicle  system  is 
introduced,  which  allows  investigating  the  behaviour  of  the  real 
components  in  a  virtual  environment.  The  main  components  of 
the  fuel  cell  HIL  system,  their  properties  will  be  described  in 
detail,  are  the  command  station,  target  node  PC  with  real  time 
operating  system,  signal-conditioning  unit,  device  under  test 
and,  for  presentation  purposes,  a  separate  human  machine  inter¬ 
face  consisting  of  a  steering-wheel  and  a  break  and  a  throttle 
pedal. 

3.1.1.  Fuel  cell  vehicle  model 

CarSim,  the  fuel  cell  vehicle  model,  provides  the  complete 
vehicle  dynamics.  The  power  unit  consist  of  an  electric  motor,  an 
inverter  model  and  a  fuel  cell  model.  Fig.  12  gives  an  overview 
of  the  model  and  the  interactions  between  different  model  parts. 

Via  the  throttle  the  required  torque  is  predetermined.  The 
power  demand,  which  is  needed  for  the  required  vehicle  speed  is 
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calculated  from  the  vehicle  model,  the  aerodynamic  and  friction 
forces.  This  power  demand  is  provided  to  the  fuel  cell  model 
where  the  actual  power  according  to  the  prevailing  operational 
conditions  is  calculated  and  provided  via  the  inverter  model 
to  the  model  of  the  asynchronous  machine.  The  propulsion 
torque  is  transmitted  via  one-step  gear  to  the  CarSim  vehicle 
model.  Below  the  different  parts  of  the  model  are  described  in 
detail. 

3.1.2.  CarSim  vehicle  dynamics  model 

CarSim  is  a  commercial  program  for  vehicle  dynamics.  It 
provides  different  precompiled  vehicle  models,  for  example 
dynamic  link  libraries,  which  can  be  easily  included  in  the 
Simulink  models.  It  also  includes  a  scientific  plotter  and  an 
animator  for  displaying  the  results.  The  animator  is  able  to  work 
in  real-time,  so  it  is  possible  to  build  up  driving  simulators.  2D 
and  3D  tables,  unlimited  in  length,  parameterise  the  vehicle 
models. 


The  CarSim  model  includes: 

•  A  wide  range  of  controllers  for  braking,  steering  (for  example 
a  driver  model),  gear  shifting  and  for  event  handling. 

•  Wind  and  aerodynamic  effects  (six  aerodynamic  forces  and 
moments  can  be  applied  to  the  sprung  mass). 

•  3D  road  geometry  and  friction. 

•  A  full  nonlinear  kinematical  behaviour  of  the  suspension. 

•  Different  tire  models  like  Pacejka  5.2,  Shear. 

•  A  power  train  model  for  front-,  rear-  and  four-wheel 
drive. 

The  model  consists  of  10  bodies  and  27  degrees  of  freedom. 
Fig.  13  gives  an  overview  about  different  modelled  motions  of 
the  vehicle. 

The  different  bodies  are  one  sprung  body,  two  unsprung  bod¬ 
ies,  rigid  axle  stroke,  rigid  axle  rotation,  four  spinning  wheels 
and  engine  crankshaft. 


Fig.  13.  Degrees  of  freedom  of  the  CarSim  vehicle  model. 
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Fig.  14.  Characteristic  map  of  asynchronous  motor. 


If  the  vehicle  model  is  used  in  Simulink  one  can  export 
approximately  530  variables  and  import  about  150  to  expand 
the  vehicle  models  with  controllers  or  vehicle  parts.  Especially 
for  the  HIL  application  CarSim  RT  includes  models  to  be  imple¬ 
mented  into  specific  real-time  platforms,  like  dSpace,  RT-Lab, 
ADI  and  more.  A  more  in  depth  description  would  be  beyond 
the  scope  of  this  article  due  to  the  complexity  of  the  model.  An 
extensive  description  can  be  found  in  literature  [5,6]. 


3.1.3.  Motor  and  inverter  model 

By  default  CarSim  offers  the  possibility  to  choose  between 
different  internal  combustion  engines  (ICE),  which  are  rep¬ 
resented  by  their  characteristic  map.  In  order  to  use  CarSim 
together  with  a  fuel  cell  as  a  power  source,  the  ICE  is  replaced 
by  the  model  of  an  electrical  motor  and  inverter. 

The  model  of  an  asynchronous  motor  combined  with  a  stan¬ 
dardised  inverter  is  employed  because  of  the  high  development 
stage.  Characteristic  maps  as  shown  in  Fig.  14  describe  the 
behaviour  of  the  motor  and  inverter. 

Fig.  15  illustrates  how  the  power  unit  reacts  on  applying  the 
throttle  pedal  by  the  driver.  The  part  marked  blue  represents 
thereby  the  motor  model. 

The  driver  applies  the  throttle,  whereas  a  final  controlling 
element  translates  the  input  to  a  claimed  torque.  Although  the 
driver  applies  the  throttle  to  reach  a  specific  speed,  the  throttle 
control  requests  a  torque.  The  driver  acts  as  a  superposed  speed 
controller.  The  torque  of  the  engine  is  obtained  from  the  angular 
speed  of  the  engine  and  the  characteristic  curve  of  the  torque, 
Fig.  14. 

The  requested  torque  together  with  the  actual  angular  speed 
of  the  electric  drive  provides  the  requested  power.  In  the  block 
“Power  balancing”  the  power  is  limited  by  the  maximum  avail¬ 
able  power  of  the  motor.  Before  entering  the  fuel  cell  model  the 
requested  power  is  corrected  for  losses  of  inverter  and  motor 
by  the  efficiency  map.  Finally,  the  fuel  cell  stack  is  loaded  with 
this  requested  power  and  drives  the  motor.  The  model  can  eas- 
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Fig.  15.  Overview  of  the  motor  and  inverter  model. 
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Fig.  17.  HIL  application  for  a  fuel  cell  vehicle. 
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ily  be  extended  if  there  is  a  need  for  more  detailed  electrical, 
mechanical  or  thermal  models  of  the  motor  and  inverter. 

3.1.4.  One -dimensional  fuel  cell  stack  model 

The  one-dimensional  fuel  cell  model  used  in  this  work  is 
described  in  more  detail  in  references  [7-9].  Different  lay¬ 
ers  of  the  fuel  cell  (Fig.  1)  with  their  different  physical  and 
electrochemical  characteristics  are  implemented  in  the  one¬ 
dimensional  model. 

The  main  equations  describing  the  dynamic  behaviour  are 
given  in  Chapter  2.  In  contrast  to  a  described  three-dimensional 
model  here  the  discretization  is  done  only  perpendicular  to 
the  species  flow,  whereas  all  the  channels  in  the  flowfield  are 
summed  up  to  one  volume. 

4.  HIL  set-up 

The  following  chapter  gives  an  overview  of  the  HIL  system 
and  the  resulting  advantages  for  the  development  process. 

The  HIL  system  consists  of  the  command  station,  target  node 
PC  with  real-time  operating  system,  signal-conditioning  unit, 
device  under  test  and  a  separate  human  machine  interface  (HMI) 
consisting  of  a  steering-wheel  and  a  break  and  a  throttle  pedal 
(Fig.  16).  The  HMI  allows  user  to  interact  directly  with  HIL 
system  and  to  “drive”  the  fuel  cell  car  in  a  3D  environment. 

The  command  station  communicates  via  ethernet  with  the 
target  PC  where  the  vehicle  model  including  inverter  and  electric 
motor  model  is  running  in  real-time.  Via  analogue  I/Os  the  target 
PC  communicates  with  the  device  under  test  (DUT).  In  this  case 
the  DUT  plays  the  role  of  a  fuel  cell  stack  and  gas  supply  units, 
which  are  also  “simulated”  and  running  in  real-time  on  a  micro 
controller  MICROGen  MPC555.  This  arrangement  is  chosen  in 
order  to  minimise  the  development  costs  of  the  HIL  system.  In 
a  real  application  the  HIL  capable  test  bench  where  the  above- 
mentioned  components  are  real  existent  could  easily  replace  the 
micro  controller. 

MAGNUM  is  a  leading  company  in  the  area  of  HIL  capable 
test  benches  for  fuel  cells  and  components,  further  details  about 
MAGNUM’s  fuel  cell  HIL  test  benches  can  be  found  in  [9,10]. 

Between  DUT  and  target  PC  a  signal-conditioning  unit 
(Genix)  is  introduced  in  order  to  adjust  the  signal  amplitudes 
in  needed  range. 

The  HMI  enables  the  user  to  drive  a  vehicle  like  a  real  car 
with  familiar  operating  elements. 

Fig.  17  shows  the  application  of  the  HIL  concept  for  the  fuel 
cell  vehicle,  In  contrast  to  the  fuel  cell  vehicle  model  (Fig.  12), 
it  can  be  recognised  that  the  fuel  cell  stack  is  not  a  part  of 
the  model  any  more.  The  fuel  cell  stack  acts  now  as  DUT  and 
exchanges  data  with  vehicle  model  via  analogue  I/Os.  The  inter¬ 
face  between  HMI  and  the  vehicle  model  is  also  presented  in 
Fig.  17. 

System  analysis  and  optimisation  combining  the  real  exist¬ 
ing  hardware  (DUT)  with  simulation  models  can  be  done  using 
the  HIL  set-up  shown  in  Fig.  17.  It  is  obvious  that,  e.g.,  the 
power  unit,  which  is  an  integral  part  of  DUT,  can  be  tested  in 
the  simulated  vehicle  environment.  The  advantage  is  that  before 
the  design  of  packaging  the  power  unit  with  its  components  can 


Fig.  18.  HIL  application  for  a  fuel  cell  vehicle. 
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Fig.  19.  Vehicle  speed. 

be  tested  and  optimised  in  a  user-friendly  environment,  e.g., 
easy  accessibility  of  components  in  a  test  bench,  etc.  Different 
hardware  configurations  can  be  tested  and  their  influence  on  the 
dynamic  behaviour  of  the  vehicle  and  also  the  fuel  consumption 
can  be  analysed. 

Another  important  application  of  the  HIL  simulation  is  the 
model  validation  in  real  time.  With  obtained  measured  data  from 
the  real  stack,  parameter  optimisation  of  the  stack  model  can  be 
performed.  Validated  real  time  fuel  cell  system  models  can  be 
used  for  model-based  controller  development  [11]. 
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Fig.  21.  Stack  voltage  and  current. 


Increased  flexibility  of  the  development  process  is  another 
great  advantage.  Consecutively,  the  different  component  mod¬ 
els  can  be  replaced  by  their  hardware  leading  to  an  optimised 
system. 

5.  HIL  simulation  results 

Several  examples  of  results  obtained  from  the  HIL  simulation 
are  presented  below.  Fig.  18  shows  the  vehicle  and  road  set-up 
used  for  the  simulation. 

In  the  following  figures  (Figs.  19-22)  results  from  a  cer¬ 
tain  driving  cycle  are  shown.  Fig.  19  shows  the  profiles  of  the 
input  for  the  speed  controller  and  calculated  speed  by  the  vehicle 
model  (Vehicle  total  CG).  In  Fig.  20  the  actual  value  of  stack 
power  is  presented.  The  response  of  the  stack  voltage,  current 
and  temperature  are  shown  in  Figs.  21  and  22. 


6.  Summary 

The  dynamical  3D  model  has  been  developed  to  describe 
the  partial  pressure  in  the  flowfield  and  GDL  dynamically.  The 
change  of  membrane  resistance  with  dynamic  load  and  the 
resulting  local  current  density  is  evaluated  within  the  model.  The 
results  for  the  dynamic  behaviour  of  the  PEMFC  show  excellent 
conformance  with  the  measurement.  They  provide  the  bases  for 
the  HIL  system  model. 

The  presented  HIL  system  used  as  development  platform 
offers  significant  advantages  compared  to  conventional  test  sys¬ 
tems.  Using  the  presented  model,  consisting  of  vehicle  dynam¬ 
ics,  electric  motor  and  inverter  and  a  fuel  cell  stack  a  wide  range 
of  system  designs  can  be  tested. 

The  introduced  HIL  set-up  allows  to  replace  certain  compo¬ 
nent  models  by  their  hardware  and  the  whole  system  behaviour 
can  be  analysed  under  real-time  conditions.  Furthermore  the 
HIL  system  offers  the  possibility  to  run  models  in  parallel  to  the 
investigated  hardware  and  to  validate  and  optimise  the  models, 
as  well  as  to  use  the  models  as  diagnostic  tools. 
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